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I.  INTRODUCTION 

The  theorist  who  plans  to  study  the  chemical  properties  of  surfaces  from 
a fundamental  electronic  point  of  view  is  presented  with  a number  of  difficult 
problems.  For  example,  if  one  focuses  on  heterogeneous  catalysis  by  metal 
particles  dispersed  on  a support^  should  the  electronic  structure  of  the 
particles  be  viewed  in  terms  of  the  band-structure  theory  of  solid-state  physics 
or  in  terms  of  the  chemical-bond  concept  of  inorganic  chemistry?  What  are  the 
most  significant  configurations  of  adsorbates  which  are  precursors  to  hetero- 
geneous reactivity?  Does  the  supporting  material  affect  chemisorption  and 
catalytic  activity?  The  theorist  who  calculates  electronic  structure  is 
accustomed  to  having  the  positions  and  configurations  of  atoms  in  a molecule 
or  solid  as  input  to  the  calculation.  Unfortunately,  detailed  structural 
information  for  heterogeneous  catalysts  is  frequently  unavailable  and  is  just 
beginning  to  be  a subject  of  experimental  investigation.  In  contrast,  the 

molecular  structures  of  transition-metal  coordination  complexes  involved  in 

2 

homogeneous  catalysis  and  enzymes  involved  in  biocatalysis  are  often  known 
to  reasonable  accuracy. 

Given  at  least  some  knowledge  of  the  catalyst  molecular  structure,  a 
number  of  new  questions  can  be  asked.  How  mobile  or  labile  are  the  adsorbates 
which  are  reaction  precursors  on  the  catalyst  surface?  Following  the  initial 
step  of  chemisorption,  what  reaction  intermediates  or  decomposition  products 
of  the  adsorbates,  if  any,  are  important?  How  do  these  intermediates  convert 
or  interact  with  other  chemisorbed  species  to  produce  the  reaction  products? 
These  questions  are  connected  with  the  nature  of  the  reaction  kinetics  and, 
from  the  theoretical  viewpoint,  are  much  more  difficult  to  answer  than  those 
underlying  simple  chemisorption. 

Keeping  these  limitations  in  mind,  one  may  take  a more  optimistic  outlook 
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and  assess  how  theory,  specifically  quantum  chemistry,  can  currently  be  used 
to  promote  the  understanding  of  catalysis  and  other  surface  chemical  properties 
at  the  electronic  level.  It  is  essential,  first  of  all,  to  investigate  the 
electronic  structures  of  polyatomic  aggregates  and  clusters  of  the  types  that 
constitute  the  active  centers  of  commercial  heterogeneous  catalysts,  along 
with  theoretical  studies  of  ideal  crystal  surfaces  and  chemisorption  thereon. 

It  Is  also  important  to  try  to  understand  the  coordination  chemical  bonding 
of  Isolated  transition-metal  complexes  of  the  type  which  are  involved  In 
homogeneous  catalysis  and  the  electronic  structures  of  the  biocatalytlcal ly 
active  centers  of  enzymes.  As  will  be  emphasized  in  this  paper,  there  are 
striking  analogies  between  catalytical ly  active  sites  on  transition-metal 

p 

surfaces  and  coordinatively  unsaturated  transition-metal  complexes  (Section 

111),  as  well  as  analogies  between  the  active  centers  of  supported  metal 

catalysts  (Section  IV)  and  the  active  centers  of  certain  metal loenzymes 

(Section  V).  Such  analogies  are  probably  not  fortuitous.  They  should  be 

investigated  and  a common  basis  of  understanding  established. 

In  this  paper  we  focus  on  a number  of  important  similarities  among  the 

electronic  structures  of  heterogeneous,  homogeneous,  and  biological  catalysts, 

using  the  results  of  recent  self-consistent-field  X-alpha  scattered-wave 

3 

(SCF-Xa-SW)  molecular-orbital  calculations.  Earlier  phases  of  this  work 

4 5 

are  described  In  previous  publications.  ’ One  of  the  main  objectives  of  such 
theoretical  studies  is  the  extraction  of  catalytic  indices,  dependent  on  the 
electronic  structure  of  the  catalyst,  which  would  be  an  approximate  gauge  of 
the  relative  activity  and  selectivity  of  the  catalyst  with  respect  to  certain 
reactants.  The  establishment  of  such  indices  could  ultimately  serve  as  a 
guide  in  the  molecular  design  of  catalysts  for  specific  reactions.  Recent 
progress  in  the  development  of  such  indices  is  the  subject  of  the  following 
section. 
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II.  SPIN-ORBITAL  ELECTRONEGATIVITY  ANO  CATALYTIC  ACTIVITY 

An  outgrowth  of  this  theoretical  work  is  the  establ ishment  of  "sjain- 

orbital  electronegativity,"  defined  by  the  SCF-Xa  spin-orbital  energies  for 

representative  surface  clusters  and  coordination  complexes,  as  a reactivity 

Index  of  the  type  described  in  the  preceding  section.  The  concept  of  spin- 

orbital  electronegativity  Is  derived  from  the  fact  that  the  orbital  energy 

eigenvalues  in  the  SCF-Xa  theory  are  rigorously  equal  to  first  derivatives 

3 6 

of  the  total  energy  with  respect  to  orbital  occupation  number,  ’ i.e., 

E1lfa  • 3<EXa>/8n1  (1) 

These  quantities  should  not  be  identified  with  the  orbital  energies  defined 

In  conventional  Hartree-Fock  theory,6  namely,  as  the  differences 

e1HF  = <EHF(ni=1)>  " <EHF(ni=0)>  (2) 

between  single-determinant  total  energies  calculated  when  the  ith  orbital  is 

occupied  and  when  It  is  empty  (fixing  the  remaining  occupied  orbitals). 

The  Xa  orbital  energies  defined  in  Eq.  (1)  correspond  closely  to  the  orbital 


electronegativities 


Xi  = 8E/3ni 
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defined  by  Hinze  e_t  aj_/  as  a generalization  of  Mulliken’s  definition  of 


electronegativity 


X„  = ]{!*A) 


where  I is  the  Ionization  potential  and  A is  the  electron  affinity  of  a 
chemically  bonded  atom  in  its  valence  state. 

This  generalization,  its  relationship  to  SCF-Xa  theory,  and  ultimately 
Its  use  as  an  Index  of  reactivity  can  be  understood  better  if  one  recalls 

9 

that  electronegativity,  as  originally  defined  by  Pauling  is  a measure  of 
the  power  of  a chemically  bonded  atom  to  attract  electrons  to  itself. 

Pauling  believed  that  electronegativity  is  a virtually  constant  atomic 
property,  even  for  different  valence  states  of  the  same  element,  and  establ ished 
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a scale  of  electronegativities  for  the  elements  based  on  the  empirical  bond 
energies  of  heteronuclear  diatomic  molecules.  Despite  the  arbitrariness  of 
this  scale  and  the  uncertainties  in  the  thermochemical  data  on  which  it  is 
based,  a wide  variety  of  chemical  phenomena  have  been  reasonably  explained 
through  the  use  of  Pauling's  electronegativity  scale. 

Pauling's  concept  of  electronegativity  as  a fixed  atomic  characteristic 
Is  somewhat  more  restrictive  than  Mulliken's  definition  of  electronegativity 
In  terms  of  the  valence-state  ionization  potential  I and  electron  affinity  A, 
since  one  does  not  expect  that  I and  A or  their  average  should  be  constant 
for  different  valence  or  oxidation  states  of  a chemically  bonded  atom.  Thus 
while  Mulliken's  electronegativity  scale,  in  its  simplest  form,  can  be 
adjusted  to  agree  reasonably  well  with  Pauling's  scale,  element  by  element, 
Mulliken's  concept  is  more  satisfying  from  a theoretical  point  of  view  and 
allows,  in  principle,  for  the  dependence  of  electronegativity  on  the  chemical 
environment  of  an  atom. 

Since  I and  A are  quantities  related  respectively  to  the  removal  of  an 
electron  from  the  highest  occupied  atomic  orbital  and  the  addition  of  an 
electron  to  the  lowest  unoccupied  orbital,  it  might  be  expected  that 
Mulliken's  concept  of  electronegativity  could  be  further  generalized  to  all 
the  orbitals  of  a chemically  bonded  atom  and  indeed  to  the  molecular  orbitals 
of  an  aggregate  of  atoms.  Thus  one  is  led  to  the  concept  of  orbital 
electronegativity  as  a measure  of  the  power  of  a chemically  bonded  atom  or 
molecular  aggregate  to  attract  an  electron  to  a particular  atomic  or  molec- 
ular orbital.  The  mathematical  definition  of  orbital  electronegativity  as 
the  first  derivative  of  the  total  energy  with  respect  to  occupation  number, 
given  In  Eq.  (3)  In  the  form  suggested  by  Hinze  et  al_.  J is  consistent  with 
the  above  conceptual  definition.  Implicit  in  Eq.  (3)  are  the  two  assumptions 


Is  a continuous  and  differentiable  function  of  the  occupation  numbers. 


The  occupation  numbers  n^  and  statistical  total  energy  E^  defined  in 
the  Xa  density* functional  self-consistent-field  theory  fulfill  both  of  the 
above  conditions,  so  that  one  can  uniquely  identify  the  SCF-Xa  electronic 

1 

energy  eigenvalues  of  an  atom,  molecule,  or  cluster,  as  given  in 

expression  (1),  with  the  orbital  electronegativities  defined  in  Eq.  (3). 

In  the  limit  where  the  total  energy  is  a quadratic  (parabolic)  function  of 

the  occupation  number,  the  orbital  electronegativity  reduces  exactly  to 

Mulliken's  definition  of  electronegativity  given  in  Eq.  (4).^  This  follows 

from  a simple  geometric  theorem  which  relates  the  slope  of  the  chord  of  a 

parabola  to  the  slope  of  the  parabola  at  its  midpoint.  The  same  type  of 

argument  applied  to  SCF-Xa  orbitals  leads  to  Slater's  transition-state 

concept,6  whereby  one  determines  1 or  A for  an  atom,  molecule,  or  cluster 

by  subtracting  or  adding  one-half  a unit  of  valence  orbital  electronic  charge 

and  then  calculating  self-consistently  the  energy  of  the  relaxed  orbital. 

While  these  relaxed  transition-state  energies  can  be  individually  identified 

with  the  corresponding  orbital  ionization  potentials  or  electron  affinities, 

the  unrelaxed  ground-state  SCF-Xa  orbital  energies  define  a set  of 

orbital  electronegativities.  Thus  the  relative  positions  of  the  SCF-\o 

electronic  energy  levels  for  a system  of  Interacting  or  reacting  atoms, 

molecules,  or  clusters  are  a measure  of  the  orbital  electronegativity  and 

chemical -potential  differences  between  the  various  reactants. 

In  those  systems  where  magnetic  spin  polarization  is  important,  one  may 

use  the  spin-unrestricted  version  of  the  SCF-Xa  method  to  calculate  different 

orbitals  for  different  spins,  leading  to  spin-polarized  energy  levels  i . ^ t- 
3 6 

and  * If  these  spin-dependent  orbital  energies  are  identified 

with  orbital  electronegativities,  then  one  is  automatically  led  to  the  concept 
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of  s pin -orb i tal  electronegativity  as  a measure  of  the  power  of  an  atom  or 
molecular  aggregate  to  attract  an  electron  to  a particular  atomic  or  molec- 
ular spin  orbital.  For  example,  the  spin  dependence  of  orbital  electro- 
negativity is  central  to  understanding  the  catalytic  activity  and  surface 
reactivity  of  iron  (see  Section  TV). 

To  understand  the  relationship  between  spin-orbital  electronegativity, 

as  represented  by  the  SCF-Xa  spin-orbital  energy  levels  of  representative 

surface  clusters  and  reactants,  and  catalytic  activity,  it  is  helpful  to 

recall  the  following  concept  originally  introduced  by  Fukui.10  For  a 

concerted  chemical  reaction  to  occur  with  reasonable  activation  energy, 

electrons  must  be  able  to  flow  between  the  reactants  from  occupied  orbitals 

into  unoccupied  orbitals  with  which  they  have  net  positive  overlap,  as  the 

reactants  move  along  the  reaction  coordinate.  Overlap  and  electron  flow 

will  be  ensured  if  the  pertinent  reactant  orbitals  have  the  following 

characteristics:  (a)  the  same  symmetry  (i.e.,  orbital  symmetry  conservation 

as  originally  proposed  by  Woodward  and  Hoffmann^1),  and  (b)  equal  or  nearly 

equal  orbital  electronegativities.  In  the  limit  where  electron  flow  between 

reactants  is  simply  from  the  highest  occupied  molecular  orbital  (HOMO)  to 

the  lowest  unoccupied  molecular  orbital  (LUMO),  condition  (b)  is  equivalent 

to  the  requirement  that  the  energy  gap  between  HOMO  and  LUMO  be  as  small 
12 

as  possible. 

Where  direct  electron  flow  between  reactants  is  forbidden  by  orbital 
symmetry  restrictions  or  unfavorable  orbital  electronegativity  differences 
(implying  a large  activation  energy),  a surface  can  heterogeneously  catalyze 
the  reaction  by  providing  a pathway  for  such  electron  flow,  e.g.,  through 
chemisorption  via  spatially  directed  or  hybridized  d-orbitals  in  the  case 
of  a transition-metal  surface.  Similar  arguments  are  applicable  to  isolated 
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transition -metal  coordination  complexes  'which  homogeneously  catalyze  electron 
flow  between  reactants  through  bonding  and  exchange  of  ligands. 

III.  COORDI NATIVELY  UNSATURATED  TRANSITION-METAL  COMPLEXES 

AS  ANALOGUES  OF  CATALYTICALLY  ACTIVE  SITES  ON  TRANS  IT I0N-MFTAL  SURFACES 

It  Is  well  known,  for  example,  that  certain  "coordinati vely  unsaturated" 

transition -metal  complexes  in  solution  can  homogeneously  catalyze  chemical 

reactions  ,2,1^h11e  It  has  long  been  suspected  that  low-coordination  sites 

on  transition-metal  surfaces  and  supported  transition-metal  clusters  are 

14  is 

centers  of  heterogeneous  catalytic  activity.  ’ ' In  this  section,  it  will 
be  shown  that  the  SCF-Xa  electronic  structure  of  such  complexes,  in  conjunc- 
tion with  the  concept  of  orbital  electronegativity,  is  consistent  with  their 
reactivity  and  is  suggestive  of  how  low-coordination  sites  on  transition- 
metal  surfaces  can  act  as  centers  of  catalytic  activity.  The  dissociation 
and  reactivity  of  is  considered  as  an  illustrative  example. 

As  a working  model,  we  consider  a Group-VIll  transition-metal  atom  (M) 
dihedrally  coordinated  by  ligands  (L),  yielding  the  coordinati vely  unsaturated 
L2M  complex  Illustrated  at  the  top  of  Fig.  1.  This  model  has  the  advantage 
that  it  can  realistically  represent  transition-metal  complexes  of  the  type 

(e.g.,  M ® Pt,  Ir,  Rh;  L = Ph^P  * triphenylphosphine)  which  dissociatively 

1 ^ 

bind  and  homogeneously  catalyze  reactions  of  Hj>,  ' and  it  can  simulate  low- 
coordination  sites  (e.g.,  "corner  atoms")  of  faceted  transition-metal  clusters 
or  stepped  transition-metal  surfaces  which  dissociati vely  chemisorb  and 
heterogeneously  catalyze  reactions  of  In  the  latter  systems,  the 

ligand  (L)  is  also  a metal  atom,  either  of  the  same  species  as  the  transition 
metal  (M),  or  of  a different  species  in  the  case  of  an  alloy  surface  or 
bimetallic  cluster. 

Molecular-orbital  calculations  have  been  carried  out  for  L-,M  and  L -.MM , 

C C 

complexes  by  the  SCF-Xa-SW  method  as  a function  of  metal  species  (M  = Pt,  1 r ) , 
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ligand  species  (L  * phosphine,  Pt),  and  molecular  geometry.  The  resulting 
orbital  energies  for  M = Pt,  L = phosphine,  and  geometry  characteristic  of 
the  platinum-phosphine  complexes  described  in  Ref.  13  are  shown  in  Fig.  1. 

Also  shown,  for  comparison,  are  the  SCF-Xa  orbital  energies  for  the  isolated 
metal,  ligand,  and  hydrogen  molecule  at  the  free-molecule  internuclear 

P O 

distance  0.74  A (F^)  and  internuclear  distance  2.8  A (F^*)  characteristic 
of  the  partially  dissociated  ("di hydride")  configuration  of  F^  in  the  l^MF^ 
complex. 

The  SCF-Xa  orbital  energy  eigenvalues  shown  in  Fig.  1 can  be  rigorously 
identified  with  orbital  electronegativities  which  are  a measure  of  the 
relative  average  electron  donor-acceptor  character  of  the  individual  orbitals, 
as  described  in  the  preceding  section.  Thus  the  fact  that  the  isolated  ligand 
energy  level,  which  corresponds  to  a phosphine  "lone-pair"  orbital,  nearly 
coincides  with  the  d-orbital  energy  level  of  the  isolated  Pt  atom  (neglecting 
relativistic  shifts)  implies  a predominantly  covalent  L-Pt(5d)  interaction 
similar  to  that  expected  for  a direct  Pt(5d)-Pt(5d)  interaction.  In  this 
respect,  the  effect  of  coordinatively  unsaturated  phosphine  ligands  on  the 
electronic  structure  of  a platinum  atom  is  expected  to  be  similar  to  that  of 
embedding  a Pt  atom  in  a low-coordination  Pt  environment,  such  as  that  provided 
by  a surface  or  cluster. 

The  ligand-metal  interaction  in  the  complex  leads  to  the  bonding 
orbital  energies  labeled  L-M(d^2)  and  L-M(d22)  in  Fig.  1,  and  to  the  anti- 
bonding  orbital  energies  labeled  Mfd^J-L*,  M(d^z)-L*,  and  M(s)-L*,  of 
which  M(dy2)-L*  is  the  highest  occupied  energy  level  in  the  ground  state  of 
the  complex.  A simple  interpretation  of  the  position  of  the  latter  energy 
level  is  that  the  strong  ligand-field  repulsion  of  the  metal  d-orbital  pointed 
on  the  ligand  directions  (the  d^2  orbital  for  the  chosen  coordinate  system) 
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raises  the  energy  level  of  this  orbital,  reduces  the  corresponding  orbital 
electronegativity,  and  mixes  in  significant  antibonding  ligand  character. 

The  dz2  orbital  is  also  subject  to  some  antibonding  ligand-field  repulsion, 
whereas  the  dx2_y2»  dxz*  anc*  °*Xy  orbitals  remain  essentially  nonbonding. 

When  platinum  atoms  are  substituted  for  the  phosphine  ligands,  the  electronic 
structure  reduces  to  the  manifold  of  bonding,  nonbonding,  and  antibonding 
d-orbital  energy  levels  (the  "d-band")  characteristic  of  a small  platinum 
cluster.^  In  this  case,  the  M(dyZ)-L*  (L  = M)  antibonding  orbital  may  be 
interpreted  as  the  analogue  of  a localized  "surface  state"  which  is  split  off 
from  the  top  of  the  d-band. 

The  most  important  result  of  the  strong  ligand-metal  anti  bonding 
component  is  to  bring  the  M(dyZ)-L*  orbital,  the  highest  occupied  orbital, 
closer  in  energy  and  electronegativity  (as  compared  with  the  isolated  Pt 
atom)  to  the  empty  antibonding  ou  orbital  of  the  molecule.  This  facilitates 
overlap  and  electron  flow  between  the  M(dyz)-L*  and  ou  orbitals,  which  are 
syirmetry  conserving,11  thereby  promoting  dissociation  of  The  partially 
di  ssociated  molecule  (H2*),  characterized  by  and  ou  orbital  energies 
approaching  the  SCF-Xa  Is  orbital  energy  of  a free  hydrogen  atom  (see  Fig.  1), 
can  bind  or  "chemisorb"  in  a dihydride  configuration  to  the  coordinatively 
unsaturated  metal  site.  This  is  revealed  by  the  L^MI^  molecular-orbital 
energies  shown  in  Fig.  1 and  the  corresponding  orbital  wavefunction  contour 
maps  shown  in  Fig.  2.  The  2b2  orbital,  for  example,  results  from  overlap 
and  electron  flow  between  the  M(d^z)-L*  orbital  and  the  ou  orbital.  The 
dihydride  configuration  is  further  stabilized  by  the  "butterfly-like"  la-j 
and  2a-j  orbitals  shown  in  Fig.  2,  formed  from  the  overlap  of  the  equatorial 
parts  of  the  L-M(dz2)  and  M(dz?)-L*  orbitals  with  the  H(ls)  (or  H?*  og) 
orbitals.  Note  that  the  M(dz2 ) lobe  pointed  along  the  z-diroct,ion  acts  as 


a repulsive  barrier  which  helps  to  keep  the  H atoms  apart.  These  dihydride 

bonding  orbitals  are  offset  somewhat  by  the  4a^  and  3a^  orbitals  resulting 
respectively  from  the  antibonding  Interaction  of  the  L-Mfd.^)  and  M( dx2_y2 ) 

orbitals  with  the  Hg*  Og  orbital,  as  is  evident  in  the  4a-j  orbital  contour 

map  shown  In  Fig.  2.  There  is  negligible  contribution  of  the  M(s)  orbital 

component  In  the  binding  of  hydrogen  to  these  platinum  and  Iridium  complexes. 

This  is  consistent  with  the  finding,  based  on  SCF-Xn  cluster  calculations^ 

and  photoemission  studies,^  that  the  metal  d-orbitals  are  almost  exclusively 

responsible  for  the  chemisorption  of  hydrogen  on  second-  and  third-row 

transition  metals  such  as  palladium  and  platinum,  whereas  significant  metal 

s,d-hybr1d1zation  (with  the  s-orbital  component  dominant)  Is  Involved  in 

hydrogen  chemisorption  on  first-row  transition  metals  such  as  nickel.  Since 

the  deuterium  molecule  (D2)  is  orbital ly  Identical  to  the  hydrogen  molecule 

(H2),  all  the  results  described  above  for  the  dissociation  of  Hg  at  a low- 

coordination  transition-metal  site  apply  equally  well  for  the  dissociation  of 

D2  at  such  a site. 

The  above-described  electronic  structure  of  the  (or  L2MD2) 

coordination  complex  leads  to  possible  explanations  of  the  observed  homo- 
geneous and  heterogeneous  catalytic  reactivity  of  H2  (or  Dg ) . For  example, 
the  near  cancellation  of  the  contributions  of  the  bonding  (la-j,  2a-|)  orbitals 
and  antibonding  (3a-|,  4a^ ) orbitals  to  metal-hydrogen  bond  strength,  leaving 
the  dissociative  2b2  bonding  orbital  dominant,  explains  the  relatively  weak, 

reversible  binding  of  H2  (or  D2)  to  such  complexes  and  their  ability  to 

1 3 

activate  H2-D2  exchange.  Since  such  a complex  is  also  a good  model  for  H2 
(or  Dg)  dissociation  at  the  corner  atoms  of  a platinum  surface  step,  the 
results  suggest  why  atomic  steps  on  platinum  surfaces  can  be  active  sites  for 
dissociating  H2  and  D2  and  in  activating  H2-D2  exchange. 

The  electronic  structure  of  the  complex  also  suggests  a possible 
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reaction  path  for  the  hydrogenation  of  unsaturated  hydrocarbons  at  low- 
coordination  transition-metal  sites.  The  4a-j  orbital,  which  defines  the 
Fermi  energy  of  the  site.  Is  closely  matched  In  symmetry,  energy,  and  electro- 
negativity to  the  it  orbitals  of  hydrocarbons  such  as  acetylene  (C2H2)  and 
ethylene  (C2H^).  When  the  4a^  orbital,  which  Is  an  antibonding  mixture  of 
L-M(d22)  and  H2*  Og  orbital  character,  is  only  partially  occupied  (as  is 
the  case  for  M * Ir,  Rh),  It  offers  a pathway  for  electron  flow  from  a C2H2 
(or  CgH^)  it  orbital  to  the  dlssociati vely  "chemisorbed"  hydrogen.  Electron 
flow  directly  between  C2H2  (or  C2H^)  and  H2  in  the  gas  phase  via  the  filled 
it  and  Og  orbitals  Is  forbidden  by  the  Pauli  exclusion  principle,  whereas 
electron  flow  directly  between  the  it  orbital  and  empty  ou  orbital  is  forbidden 
by  orbital  symmetry. ^ Because  the  4a^  orbital  of  L2MH2  is  antibonding  between 
the  L2M  site  and  H2,  while  bonding  between  H2  and  C2H2  (or  C^),  the  net 
result  of  electron  flow  between  a n orbital  and  the  4a^  orbital  is  the 
breaking  of  a C-C  tt  bond,  the  formation  of  two  new  C-H  bonds,  and  the  expulsion 
of  the  hydrogenated  species  C2H^  (or  C2Hg),  as  suggested  by  the  reaction  path 
shown  In  Fig.  3.  Also  shown  is  a contour  map  for  the  4a^  orbital  of  the 
L2MH2C2H2  reaction  Intermediate  (the  third  step  of  the  proposed  reaction  path) 
formed  as  a result  of  the  Interaction  of  acetylene  with  the  L2MH2  complex. 

The  Incipient  formation  of  C-H  bonds  via  the  overlap  of  the  C-C  n orbital 
with  the  antibonding  me tal-dl hydride  orbital  and  the  resulting  ethylene-1  ike 
configuration  are  clearly  visible  in  this  map.  It  Is  important  to  note  that 
the  concerted  reaction  path  Indicated  in  Fig.  3 is  not  the  conventional  one 
for  hydrogenation  on  Ideal  transition-metal  surfaces,  where  it  is  usually 
assumed  that  chemisorption  of  acetylene  or  ethylene  on  one  or  two  metal 
sites  Is  the  precursor  to  combining  with  hydrogen  chemisorbed  on  neighboring 
sites.  Nonconccrted  reaction  paths  in  which  both  reactants  are  coordinated 
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to  the  same  metal  site  are  also  possible  and  indeed  have  been  argued  to  be 

1 s 

favored  kinetically  in  certain  homogeneous  reactions.  Alternative  reaction 
paths  at  low-coordination  transition-metal  sites  are  currently  under  investi- 
gation in  conjunction  with  theoretical  studies  of  the  reactivity  of 
IrCl(C0)(Ph3P)2  (Vaska's  complex).^ 

In  this  section,  we  have  attempted  to  show  that  a detailed  theoretical 
study  of  the  electronic  structure  of  well  characterized  coordinatively 
unsaturated  transition-metal  complexes  and  their  interactions  with  H2  can 
not  only  lead  to  an  understanding  of  their  homogeneous  reactivity  but  can 
also  serve  as  a model  for  the  dissociative  chemisorption  and  heterogeneous 
reactivity  of  H2  on  low-coordination  transition-metal  surface  sites,  where 
definitive  structural  information  is  lacking.  There  are  many  other  useful 

analogies  which  can  be  made  between  molecular  transition-metal  coordination 

20  21 

complexes  and  surface-adsorbate  interactions.  * 

In  concluding  this  section,  it  is  important  to  compare  the  present 
theoretical  approach  to  surface  reactivity,  based  on  molecular-orbital 
indices,  with  other  theoretical  approaches  to  this  problem  and  to  chemical 
reactivity  in  general.  In  the  applications  of  traditional  methods  of  quantum 
chemistry  (e.g.,  Hartree-Fock,  Configuration-Interaction,  Valence-Bond,  etc.) 
to  reaction  kinetics  and  thermochemistry,  one  usually  focuses  directly  on 
the  calculation  of  total  energies,  total -energy  differences,  and  "potential 
surfaces"  for  the  reactants.  This  is  generally  a computationally  difficult 
and  costly  process  to  carry  out  over  the  various  possible  reaction  paths, 
even  for  the  simplest  reactions.  To  appreciate  the  magnitude  of  this 
problem,  one  need  only  consider  the  recent  status  of  the  quantitative  first- 
principles  determination  of  the  kinetics  of  one  of  the  simplest  gas-phase 

22 

chemical  reactions,  namely,  hydrogen-deuterium  exchange,  H2  + 0 ■*  HD  + H. ' 
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This  situation  hardly  makes  one  confident  In  the  efficacy  of  quantum  theory 
to  "predict"  the  path  or  kinetics  of  surface-activated  reactions  on  the  basis 
of  potential-surface  computations.  Even  if  one  carries  out  total  energy 
calculations  for  only  a few  representative  molecular  configurations,  rather 
than  for  the  entire  potential  surface,  there  is  still  the  uncertainty  associated 
with  the  direct  subtraction  of  two  total  energies,  which  are  usually  large 
numbers,  to  obtain  a relatively  small  energy  difference  of  chemical  signifi- 

1 

cance. 

The  theoretical  approach  described  in  this  paper,  while  not  a substitute 
for  ab  initio  potential -surface  calculations,  circumvents  many  of  the  diffi- 
culties associated  with  the  latter  approach.  By  placing  emphasis  on  the 
determination  of  molecular-orbital  indices  of  reactivity,  rather  than  total 
energies,  for  realistic  transition-metal  coordination  complexes  and  clusters 
simulating  local  bulk  and  surface  configurations,  one  retains  the  molecular- 
orbital  picture  which  chemists  have  traditionally  used  and  makes  contact 
with  the  band-structure  concept  of  solid-state  and  surface  physics. 

IV.  ELECTRONIC  STRUCTURE  AND  CATALYTIC  ACTIVITY  OF 
TRANSITION-METAL  CLUSTERS 

There  is  much  current  interest  in  the  electronic  structures  of  small 

transition-  and  noble-metal  aggregates  which  constitute  the  active  centers 

1 15  21  ° 

of  heterogeneous  catalysts.  ’ ’ ' Metal  aggregates  less  than  10  A in  size, 

referred  to  as  "clusters"  in  order  to  distinguish  them  from  larger  particles 

24 

or  crystallites,  are  of  special  importance  because  their  electronic 

structures,  catalytic  properties,  and  interactions  with  supporting  environments 

can,  in  principle,  deviate  from  those  characteristic  of  the  bulk  metals. 

4 16 

In  recent  publications,  ’ the  electronic  structures  of  small  copper, 
nickel,  palladium,  and  platinum  clusters,  as  calculated  by  the  self-consistent- 
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field  X-alpha  scattered-wave  (SCF-Xa-SW)  molecular-orbital  method,  have  been 
described  and  compared  with  the  results  of  similar  calculations  by  semiempirical 
molecular-orbital  methods.  One  is  at  first  struck  by  the  systematic  similari- 
ties, rather  than  the  differences,  between  the  SCF-Xa  electronic  structures 
of  the  clusters  and  those  of  the  corresponding  bulk  metals.  For  example,  the 
manifolds  of  orbital  energy  eigenvalues  for  13-atom  clusters  having  the  cubo- 
octahedral  nearest-neighbor  coordination  characteristic  of  the  fee  lattice 
exhibit  all  the  principal  features  of  the  bulk  band  structures,  e.g.,  overlap 
of  the  "d  band"  by  the  "s,p  band,"  a sharp  peak  in  the  density  of  states  around 
the  Fermi  energy  for  the  Ni13>  Pd13,  and  Pt13  clusters,  increasing  band  width 

4 

through  this  series,  and  magnetic  spin  polarization  in  the  case  of  Ni'13. 

More  recent  SCF-Xa  studies  of  4-  and  6-atom  nickel,  palladium,  and  platinum 
clusters,  including  relativistic  effects,  indicate  that  even  these  small 
clusters  exhibit  most  of  the  qualitative  features  of  the  crystalline  band 
structure  and  can  be  utilized  as  a basis  for  understanding  the  nature  of  the 
interaction  of  hydrogen  with  these  metals.^  These  predictions  have  been 
confirmed  by  spectroscopic  studies  of  small  noble-  and  transition-metal 

or  nc 

clusters  isolated  in  inert  matrices.  ’ 

27 

Concurrent  with  these  theoretical  studies,  Messmer  and  Salahub  have 
carried  out  SCF-Xa-SW  calculations  for  aluminum  clusters  (and  oxygen  chemi- 
sorption thereon)  containing  up  to  43  atoms.  The  results  indicate  that  25 
atoms  are  sufficient  to  yield  the  detailed  structure  and  92%  of  the  band  width 
Of  the  electronic  density  of  states  of  crystalline  aluminum,  as  measured  by 
X-ray  and  photoelectron  emission  spectra.  Increasing  the  cluster  size  to 
43  atoms  yields  99%  of  the  crystalline  band  width.  The  studies  are  perhaps 
the  most  definitive  test  of  the  "real-space"  cluster  representation  of  solids 
and  surfaces  made  thus  far,  since  aluminum  is  traditionally  viewed  as  a highly 
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delocalized  nearly- free-electron  metal  for  which  the  conventional  "k-space" 
representation  of  solid-state  physics  is  most  appropriate. 

While  a small  transition-metal  cluster  does  indeed  exhibit  most  of  the 

characteristics  of  the  crystalline  band-structure,  the  molecular  boundary 

conditions  and  the  fact  that  most  of  the  atoms  of  the  cluster  are  effectively 

coordinatively  unsaturated  "surface"  atoms  lead  to  some  important  additional 

4 

features  of  the  electronic  structure.  These  include  the  appearance  of 
localized  orbitals  split  off  in  energy  from  the  top  and  bottom  of  the  cluster 
d-orbltal  manifold,  which  can  be  interpreted  as  the  cluster  analogues  of  the 
"surface  states"  split  off  from  the  bulk  d band  for  an  extended  crystalline 
surface.  Also  associated  with  the  cluster  boundary  conditions  is  a buildup 
of  electronic  charge  density  in  the  cluster  interior,  compensated  for  by  a 
depletion  of  electronic  charge  at  the  cluster  periphery.  Since  a faceted 
polyhedral  (e.g.,  cubo-octahedral  or  icosahedral)  cluster  is  the  simplest 
globular  analogue  of  a "stepped"  crystalline  surface,  the  effective  positive 
charge  on  the  atoms  at  the  cluster  periphery  is  analogous  to  the  positive 
charge  observed  at  the  catalytical ly  active  edge  atoms  of  a stepped  transition- 
metal  (e.g.,  platinum)  surface.  The  resulting  electric  field  gradient 
at  the  cluster  boundary,  coupled  with  a high  density  of  d orbitals  spatially 
directed  away  from  the  coordinatively  unsaturated  atoms  at  the  cluster  periphery 
(see  Fig.  17  of  Ref.  4),  can  promote  interaction  and  overlap  of  these  orbitals 
with  symmetry-conserving  orbitals  of  adsorbates,  catalytic  reactant  molecules, 
and  supporting  environments.  Thus  the  electronic  structures  of  small 
transition-metal  clusters  not  only  mimic  the  band  structures  of  the  corres- 
ponding crystalline  metals,  but  also  have  distinct  features  that  may  be  key  to 
understanding  the  active  centers  of  heterogeneous  catalysts  and  catalytically 
active  sites  (e.g.,  "steps")  on  otherwise  ideal  crystalline  surfaces. 

Although  nickel,  palladium,  and  platinum  clusters  of  the  type  studied 
in  Refs.  4 and  16  arc  a logical  starting  point  in  the  fundamental  investigation 
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of  catalysts,  other  Group-VIlI  transition  metals  are  of  comparable  catalytic 

versatility  and  invite  theoretical  consideration.  Among  these,  iron  is  the 

most  widely  used  commercial  catalyst  for  the  synthesis  of  ammonia  (the  Haber 
29  30 

process),  ’ is  used  commercially  as  a catalyst  in  the  Fischer-Tropsch 

31 

synthesis  of  high-molecular-weight  paraffins,  and  is  generally  one  of 
the  most  surface-active  transition  metals.  The  electronic  structures  and 
related  properties  of  iron  clusters  are  also  of  intrinsic  interest  in  the 
ways  they  compare  with  the  band  structure  and  physical  properties  of  bulk 
crystalline  iron.  Of  particular  importance  are  the  magnetic  states  of  iron 
clusters,  their  relation  to  the  ferromagnetism  of  crystalline  iron,  and  the 
influence,  if  any,  of  such  magnetism  on  the  catalytic  activity  and  surface 
reactivity  of  this  metal. 

In  its  biological  states,  iron  is  also  a center  of  biocatalytic  activity. 

For  example,  iron-porphyrin  complexes  are  the  active  centers  of  hemoglobin, 

32 

cytochrome,  and  other  hemoproteins  key  to  all  respiration  and  metabolism, 
whereas  Fe^Sg  clusters  (see  Fig.  8)  are  the  active  centers  of  ferredoxin, 
an  iron-sulfur  protein  that  is  believed  to  be  involved  in  nitrogen  fixation 
(the  biological  analogue  of  anmonia  synthesis).  The  fact  that  an  iron- 
sulfur  cluster  can  be  a center  of  biocatalytic  activity  is  intriguing,  since 
sulfur  is  well  known  to  be  a poison  of  commercial  iron  catalysts. 

This  section  is  devoted,  in  part,  to  a discussion  of  the  magnetic  states 
of  iron  clusters,  especially  as  a basis  for  understanding  the  local  bulk  and 
surface  magnetic  properties  of  crystalline  iron,  including  the  transition 
from  ferromagnetism  to  paramagnetism.  Also  discussed  is  the  possible 
relationship  between  iron  cluster  electronic  structure  and  catalytic  activity, 
utilizing  the  concept  of  spin-orbital  electronegativity  described  in 
Section  II.  Finally,  in  Section  V the  electronic  structure  of  iron-sulfur 
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clusters  in  ferredoxin  is  described  and  compared  with  the  electronic  structure 

of  elemental  Iron  clusters.  Apart  from  the  biological  significance  of  these 

studies,  the  results  show  the  importance  of  Jahn-Teller  distortions  in  small 

clusters,  give  Insight  as  to  how  sulfur  can  promote  the  biocatalytic  activity 

of  iron  despite  the  fact  that  sulfur  is  usually  an  iron  catalyst  poison,  and 

provide  a structurally  well  characterized  molecular  analogue  of  a metal- 

catalyst-support  system  including  the  effects  of  catalyst-support  interactions. 

Extensive  applications  of  the  SCF-Xa-SW  method  to  hemoglobin  have  recently 

34 

been  made  by  Case  and  are  reported  elsewhere. 

SCF-Xa-SW  molecular-orbital  calculations  have  been  carried  out  for  a 

35 

variety  of  iron  clusters  ranging  in  size  from  4 atoms  to  15  atoms.  In  this 
section  we  focus  on  the  results  for  9-  and  15-atom  clusters  having  the  body- 
centered-cublc  (bcc)  geometry  characteristic  of  the  local  atomic  arrangement 
In  bcc  crystalline  a-iron.  The  calculations  were  implemented  in  spin-unrestricted 
form,  l.e.,  different  orbitals  for  different  spins,  in  precisely  the  same  fashion 
as  described  for  nickel  clusters  in  Ref.  4.  The  resulting  spin-orbital  energy 
eigenvalues  for  the  Fe^  and  Fe^  clusters  are  shown  in  Figs.  4 and  5,  respec- 
tively, labeled  according  to  the  irreducible  representation  of  the  symmetry 
group.  An  arrow  points  to  the  highest  occupied  spin  orbital  and  defines  the 
"Fermi  energy"  of  each  cluster. 

From  the  qualitatively  similar  spin-polarized  electronic  structures  of 
the  Feg  and  Fe^  clusters,  one  can  extract  the  following  characteristics 
which  are  directly  comparable  with  the  band  structure  of  ferromagnetic 
crystalline  a-iron  calculated  by  Tawil  and  Callaway  using  an  SCF-Xa-LCAO 
technique: 

(1)  The  cluster  electronic  structures  can  be  characterized  in  terms  of 

manifolds  of  predominantly  <M  ike  spin  orbitals  bracketed  in  energy  by 
orbitals  which  have  significant  s.p  character  (as  labeled  in  fig.  7), 
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In  analogy  to  the  overlap  of  the  d hand  by  the  s,p  band  In  bulk  crystal- 
line iron.36 

(2)  The  Fermi  level  of  each  cluster  passes  through  the  center  of  the  "spin- 
down"  (i)  d-orbital  manifold  in  a fashion  similar  to  the  intersection  of 
the  Fermi  level  with  the  center  of  the  "minority-spin"  d band  of  ferro- 
magnetic iron  (see  Fig.  2 of  Ref.  36). 

(3)  The  exchange  splitting  and  widths  of  the  cluster  d bands,  particularly 
those  for  Fe.^,  are  comparable  with  the  values  for  bulk  ferromagnetic 
Iron. 

(4)  The  cluster  total  density  of  states  shows  a pronounced  two-peak  structure, 
which  is  primarily  a consequence  of  the  exchange  splitting,  consistent 
with  the  density  of  states  for  ferromagnetic  iron  (see  Fig.  5 of  Ref.  36). 

(5)  Bonding  and  antibonding  cluster  spin  orbitals  of  t2g  and  e^  symmetry, 
such  as  those  shown  in  Fig.  6,  have  a spatial  character  similar  to  the 
energy  eigenstates  of  ferromagnetic  iron  near  the  bottom  and  top, 
respectively,  of  the  d bands. 

(6)  The  ?t2g+  an<*  orbitals  are  the  lowest  unoccupied  majority-spin 

eigenstates  of  the  Fe^  cluster  (see  Fig.  5)  and  are  analogous  to  the 
majority-spin  crystal  eigenstates  of  ferromagnetic  iron  lying  immediately 
above  the  Fermi  energy  (see  the  band  labeled  in  Fig.  1 of  Ref.  36). 

SCF-Xa  calculations  have  also  been  carried  out  for  a bcc  Fe^  cluster 

in  which  the  boundary  conditions  have  been  modified  to  simulate  the  "embedding" 

35 

of  the  cluster  in  an  extended  crystalline  environment.  Although  the  corres- 
pondence between  the  cluster  and  crystalline  electronic  structures  is 

somewhat  improved,  e.g.,  through  an  effective  increase  in  the  cluster  d-band 

width,  the  main  conclusions  given  above  are  essentially  unaltered. 

The  cluster  results  can  be  further  utilized  as  a basis  for  understanding 

37 

the  transition  between  collective  ferromagnetism  (or  "superparamagnetism"  ) 
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and  local  paramagnetism  with  decreasing  particle  size  or  increasing  tempera- 
ture. For  example: 

(1)  The  average  paramagnetic  magneton  numbers  per  atom  for  the  Feg  and 

Fej5  clusters  are  2.9  and  2.5,  respectively,  approaching  the  2.2  value 

for  ferromagnetic  crystalline  a-iron,  whereas  SCF-Xa  results  for  smaller 

Iron  clusters,  such  as  Feg  and  Fe^,  suggest  increasing  magnetic  moment 

approaching  the  atomic  limit.  The  increasing  magneton  number  with 

decreasing  iron  cluster  size  is  supported  experimentally  by  Mossbauer 

38 

studies  of  small  iron  aggregates  isolated  in  inert  matrices  and  may 
be  contrasted  with  the  decreasing  magneton  number  with  decreasing 
cluster  size  observed  theoretically  and  experimentally  for  nickel  (see 

Ref.  4). 

(2)  The  partial-wave  decomposition  of  the  Feg  and  Fc15  cluster  spin  orbitals 

indicates  that  the  contribution  of  the  4s-like  components  to  spin 

polarization,  although  relatively  small  in  magnitude,  is  opposite  in 

direction  (antiparallel  to)  that  of  the  dominant  3d- like  components, 

in  good  agreement  with  the  analysis  of  magnetic  form  factors  in 

39 

polarized  neutron  scattering  data  for  crystalline  iron. 

(3)  Somewhat  more  spin  density  in  the  Fcg  and  Fe^g  clusters  is  concentrated 
in  the  e^  orbitals  than  in  the  orbitals  (see  Fig.  6),  providing  a 
local  model  for  the  concentration  of  spin  density  along  the  [100] 
direction  (the  direction  of  easy  magnetization)  in  crystalline  a-iron. 
as  deduced  from  neutron-diffraction  measurements. 

(4)  Neutron-diffraction  studies  of  crystalline  iron  at  temperatures  greater 

41  -43 

than  the  Curie  temperature  Tc  suggest  that  short-range  ordering 

O 

Of  spins  in  the  form  of  "spin  clusters"  ('v  10  A in  size)  persists  well 
Into  the  paramagnetic  region,  and  that  the  paramagnetic  magneton  number 
is  approximately  25%  greater  than  the  value  for  ferromagnetic  iron 
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(T  < T 1.  It  is  obvious  therefore  that  the  electronic  structure 
c 

of  crystalline  iron  in  the  paramagnetic  region  is  more  appropriately 
represented  statistically  by  a local  spin-polarized  cluster  model, 
than  by  the  extended  £-space  band-structure  description,  since  the 
cluster  electronic  structure  is  not  dependent  on  the  assumption  of 
long-range  crystalline  ordering  of  spins.  Utilizing  the  bcc  Fe^5 
cluster  as  a prototype  model  for  a "spin  cluster,"  we  see  that  the 
virtual  majority-spin  orbitals,  7t2gt  and  6egt  (see  ^9-  5)*  in  the 
ground  state  (T  = 0°K)  are  within  a fraction  of  an  electron  volt  of 
the  Fermi  energy,  in  analogy  to  the  lowest  unoccupied  majority-spin 
eigenstates  of  ferromagnetic  crystalline  iron  [see  item  (6)  described 
above].  Thermally  induced  "spin-flip"  electronic  excitations  from  the 
Fermi  level  (which  intersects  the  minority-spin  manifold)  to  the  7t2gt 
and  6egt  orbitals  for  T > Tc  results  in  an  effective  30%  increase  of 
cluster  magnetic  moment  and  magneton  number,  thereby  providing  a good 
model  for  the  observed  increase  in  magnetic  moment  of  crystalline  iron 
in  the  paramagnetic  region.  Such  excitations  only  slightly  alter  the 
average  exchange  splitting  and  total  density  of  states  of  the  cluster, 
consistent  with  the  apparent  persistence  of  the  local  density  of  states 
for  crystalline  iron  at  T > Tc  (as  measured,  for  example,  by  photo- 
emission), an  experimental  observation  which  heretofore  has  been 
inexplicable  within  the  framework  of  conventional  ferromagnetic  energy- 
band  theory. 

We  can  benefit  from  the  large  effective  surfaces  presented  by  the  Feg 
and  Fe^5  clusters  to  discuss  the  catalytic  activity,  selectivity,  and  general 
surface  reactivity  of  iron  aggregates,  as  has  previously  been  done  for  nickel, 
palladium,  and  platinum  in  Refs.  4 and  16.  Among  the  manifolds  of  densely 
spaced  d-orbital  eigenstates  of  these  clusters  are  levels  which  correspond 
to  antibonding  spin  orbitals  primarily  localized  on  and  spatially  oriented 
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away  from  the  cluster  periphery  or  "surface,"  especially  in  the  general 
vicinity  of  the  Fermi  energy.  The  2t^f  orbital  mapped  in  Fig.  6 is  one 
example.  Many  of  these  spin  orbitals  have  the  proper  spatial  character 
for  symmetry-conserving^  overlap  with  the  orbitals  of  certain  reactant 
molecules  and  are  likely  to  be  present  in  abundance  on  the  surfaces  of  cata- 
lytic iron  aggregates  and  particles  of  more  general  morphology.  This  argument 
is  clarified  by  a direct  comparison  in  Fig.  7 of  the  Fe^  cluster  spin-orbital 

energies  with  the  SCF-Xa  molecular-orbital  energies  of  N^,  CO,  and  O2,  the 

29  30 

molecules  which  are  key  reactants  in  iron-catalyzed  ammonia  synthesis  ’ 

31  44 

Fischer-Tropsch  synthesis,  and  surface  oxidation. 

As  described  in  Section  II,  the  SCF-Xa  spin-orbital  energy  eigenvalues 
of  a molecule  or  cluster  define  on  a one-to-one  basis  a manifold  of  electro- 
negativities, each  of  which  is  a measure  of  the  average  electron  donor  and 
acceptor  character  of  the  corresponding  spin  orbital.  Thus  the  relative 
positions  of  the  orbital  energies  of  the  Fe^  cluster  and  reactant  molecules, 
N2,  CO,  and  O2,  as  shown  in  Fig.  7,  are  a measure  of  the  differences  in 
orbital  electronegativity  among  these  components.  An  effect  of  the  large 
exchange  splitting  in  the  iron  clusters  is  to  raise  the  minority  spin  d 
orbitals  to  higher  energies  in  comparison  with  the  non-spin-polarized  limit, 
effectively  reducing  the  orbital-electronegativity  difference  and  enhancing 
covalent  overlap  between  the  cluster  HOMO  (those  symmetry-conserving  minority- 
spin  orbitals  in  the  vicinity  of  the  Fermi  energy)  and  the  LUMO  of  molecules 
such  as  Ng,  CO,  and  O2  (see  Fig.  7).  Because  the  latter  molecular  orbitals 
are  antibonding,  overlap  and  concomitant  reduction  in  the  activation  energy 
for  effective  electron  flow  between  the  iron  HOMO  and  reactant  IUM0  should 
promote  dissociation  of  the  reactant  molecules,  a precursor  to  their  cata- 
lytic reactivity.  Dissociation  of  N^,  CO,  and  O2  does  indeed  occur  on  iron 

44-47 

surfaces  at  room  temperature  and  low  coverage. 
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Also  Included  for  comparison  in  Fig.  7 are  the  SCF-Xa-SW  orbital  energies 

of  a Pt-jj  cluster  previously  shown  to  exhibit  most  of  the  characteristics  of 

the  bulk  and  surface  electronic  structures  of  crystalline  platinum  or  small 
4 

particles  thereof.  In  contrast  to  iron,  there  is  no  exchange  splitting  of 

the  platinum  cluster  d orbitals,  corresponding  to  the  nonmagnetic  state  of 

bulk  platinum,  and  the  Intrinsically  larger  electronegativity  of  the  platinum 

atom  results  in  a higher  effective  orbital  electronegativity  of  the  platinum 

cluster.  The  inclusion  of  relativistic  corrections  to  the  electronic  structure 

of  the  latter  cluster,  as  described  in  Ref.  16,  increases  its  effective  orbital 

electronegativity  even  further.  Consequently,  the  d orbitals  around  the  Fermi 

level  of  the  platinum  cluster,  although  high  in  density,  are  poorly  matched 

in  energy  and  orbital  electronegativity  to  the  LUMO  of  and  CO.  Thus  a 

high  density  of  states  around  the  Fermi  energy  is  not  a sufficient  condition 

for  catalytic  activity.  These  theoretical  results  are  consistent  with  the 

experimental  fact  that  platinum,  which  is  an  excellent  catalyst  for  some 

reactions.  Is  inactive  in  promoting  and  CO  dissociation  and  Is  a poor 

15  31 

catalyst  for  ammonia  and  Fischer-Tropsch  syntheses.  * The  above  conclusions 

about  the  relationship  between  the  orbital  electronegativities  of  transition- 

metal  clusters  and  the  conditions  for  molecular  dissociation  thereon  are 

very  similar  to  those  for  molecular  dissociation  on  transition-metal  surfaces 

arrived  at  empirically  by  BrodSn  et  al_.46 

Dissociation  of  02»  the  rate-limiting  precursor  to  oxygen  chemisorption 

and  surface  oxidation,  occurs  rapidly  at  relatively  low  temperatures  and 

44 

pressures  on  Iron  surfaces,  suggesting  low  activation  energy.  On  platinum 

48 

surfaces  this  process  seems  to  be  less  facile,  generally  requiring  more 

stringent  conditions  for  activation,  such  as  the  presence  of  surface  "steps" 

49 

which  play  the  role  of  active  sites.  Like  the  processes  of  N2  and  CO 
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dissociation  considered  above,  the  relative  activities  of  platinum  and  Iron 

for  O2  dissociation  can  be  rationalized  In  terms  of  the  respective  orbital 

electronegativities  defined  by  the  SCF-Xa  orbital  energy  eigenvalues  shown 

In  Fig.  7.  The  Fermi  level  of  the  Fe^g  cluster,  as  well  as  that  of  other 

spin-polarized  Iron  aggregates  Including  crystalline  a-iron,  lies  somewhat 

above  the  energy  of  the  partially  occupied  orbital  of  02,  Implying  that 

O2  Is  effectively  electronegative  with  respect  to  this  metal.  The  resulting 

electrophilic  addition  of  electrons  to  the  antibonding  orbital  will  promote 

O2  dissociation.  In  contrast,  the  Pt^  Fermi  level  lies  below  the  orbital, 

indicating  that  O2  Is  formally  nucleophilic  with  respect  to  this  cluster. 

Therefore,  under  normal  circumstances,  electron  transfer  between  02  and  a 

platinum  aggregate  of  "surface"  is  In  a direction  which  should  favor  bonding 

rather  than  dissociation  of  this  molecule.  Nevertheless,  a platinum  site 

can  be  "activated"  for  Og  dissociation  by  placing  the  former  In  a low- 

coordination  environment  In  which  a strong  ligand-metal  antibonding  Interaction 

raises  the  energy  (lowers  the  electronegativity)  of  the  Pt  electron  transfer 

from  the  metal  to  dioxygen.  Such  an  environment  Is  provided  by  a 

coordlnatlvely  unsaturated  platinum  complex  of  the  type  recently  investigated 

50 

via  the  SCF-Xa-SW  method  by  Norman,  or  by  the  "edge”  atoms  on  a platinum 
51 

surface  step,  as  described  In  Section  ill. 

In  an  attempt  to  elucidate  further  the  mechanism  of  ammonia  synthesis 
35 

on  Iron  catalysts,  Yang  has  also  constructed  theoretical  models  for  the 

effects  of  "promoters"  (substances  which,  added  to  the  catalyst,  enhance  its 

activity  ),  as  well  as  models  for  the  surface  nitride  which  is  formed  upon 

47 

dissociation  of  Ng  on  Iron.  The  principal  effect  of  the  promoter  is  electron 
donation  to  Iron,  which  reduces  the  electronegativity  difference  and  activa- 
tion energy  for  electron  transfer  between  the  metal  HOMO  and  N2  LUMO  (as 
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compared  with  pure  Iron),  thereby  promoting  ^ dissociation.  The  dissociative 
chemisorption  of  nitrogen  results  in  the  formation  of  a surface  nitride  with 
an  effective  negative  charge  which,  in  a reducing  atmosphere,  should 
facilitate  protonation  of  surface  nitrogen  atoms,  the  formation  of  N-H 
bonds,  and  ultimately  ammonia  synthesis.  An  analogous  argument  applied 
to  the  dissociative  chemisorption  of  CO  on  iron  leads  one  to  conclude  that 
the  formation  of  an  iron  surface  carbide  should  facilitate  C-H  bond  formation 
and  Flscher-Tropsch  synthesis. 

Details  of  these  theoretical  studies,  theoretical  models  for  ^ 
dissociation  and  ammonia  synthesis  on  ruthenium  and  osmium  catalysts,  and 
the  nature  of  catalyst-support  interactions  are  the  subjects  of  forthcoming 
papers. 

V.  THE  ACTIVE  CENTERS  OF  IRON-SULFUR  PROTEINS 

AS  BIOLOGICAL  ANALOGUES  OF  SUPPORTED  METAL  CLUSTERS 
52 

Iron-sulfur  proteins  contain  at  least  one  iron  atom  coordinated  by 

one  or  more  sulfur  ligands.  These  proteins  participate  in  many  biologically 

Important  reactions,  including  photosynthesis,  nitrogen  fixation,  and 

53 

hydroxylation  of  steroids.  Because  of  their  involvement  in  so  many 

fundamental  biological  processes,  iron-sulfur  proteins  are  obviously  of 

Interest  to  biochemists.  They  have  also  attracted  the  attention  of  many 

inorganic  and  physical  chemists  and  physicists  because  of  their  novel 

structural  and  spectroscopic  properties,  and  the  opportunity  they  provide 

54 

for  the  application  of  a wide  range  of  experimental  techniques. 

Although  the  typical  molecular  weight  of  these  proteins  is  between 
6,000  and  20,000,  there  are  only  one  or  two  iron-sulfur  units  (each  unit 
may  have  one,  two,  or  four  iron  atoms)  in  the  proteins  that  are  responsible 
for  most  of  the  observed  physical  properties  and  biological  activity. 
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These  iron-sulfur  units  are  called  "active  sites"  and  they  may  be  considered 

as  biological  analogues  of  the  active  centers  of  Industrial  supported  metal 

catalysts.  Over  the  past  few  years,  Holm  et^  aK  have  successfully  synthesized 

small-molecular-weight  analogues  of  the  iron-containing  active  sites,  corres- 

55  56  57 

ponding  to  the  1-Fe,  2-Fe,  and  4-Fe/8-Fe  proteins,  whose  spectroscopic 

and  magnetic  properties  and  X-ray  structures  parallel  closely  those  of  the 
protein  active  sites.  These  analogues  provide  realistic  models  for  theoretical 
studies  and  the  SCF-Xa-SW  method  can  readily  be  applied  to  them.  Such  studies 
for  the  1-Fe  and  2-Fe  analogues  have  been  carried  out  by  Norman  and  co- 
workers,  and  preliminary  results  for  a 4-Fe  analogue  have  been  reported 
by  Yang  et^  al 

In  this  section,  the  electronic  structure  of  the  4-Fe  protein  analogue 
[Fe^S^SCH-j)^]  is  discussed.  (S*  is  referred  to  as  "inorganic"  sulfur  in 
biochemistry  as  distinct  from  "organic"  sulfui  S.)  This  is  the  simplest 
4-iron-sulfur  protein  analogue  that  has  been  synthesized  by  Holm  and  co- 
workers, ^ and  It  Is  equivalent  to  the  active  site  of  the  reduced  high- 
potential  Iron  protein  (HPred)  and  the  active  site  of  the  oxidized  ferredoxin 
(pdox).  I.e.,  [Fc4SJ|(SCH3)4]?"  h HPrcd  = FdQX.  (The  high-potential  iron 
protein  can  be  isolated  from  the  photosynthetic  bacterium  "chromatium 

vinosum,"  and  the  ferredoxin  can  be  obtained  from  bacterium  "peptococcus 
62 

aerogenes."  ) 

A fundamental  understanding  of  the  4-iron-sulfur  proteins  may  be  key 

to  understanding  the  mechanism  of  nitrogen  fixation,  the  biological  equivalent 

33 

of  ammonia  synthesis.  Nitrogenase,  which  is  the  enzyme  complex  capable 
of  reducing  dinitrogen  to  ammonia,  consists  of  one  iron  protein  and  one 
molybdenum- Iron  protein.  It  is  believed  that  the  Mo-Fe  protein  functions  to 
bind  nitrogen,  and  the  Fe  protein  or  an  external  ferredoxin  (containing 
4-Fe-S  units)  supplies  low  potential  electrons  to  effect  the  reduction. 
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More  recently,  it  has  been  reported  J that  the  Mo-Fe  protein  also  contains 

4-Fe-S  units  which  are  very  similar  to  those  found  in  the  HP  proteins. 

From  the  point  of  view  of  catalytic  chemistry,  it  is  useful  to  study 

these  4-iron-sulfur  proteins  and  their  analogues,  not  only  because  it  may 

possibly  lead  to  a better  understanding  of  the  interaction  mechanism  between 

nitrogen  and  the  corresponding  active  centers  in  a commercial  catalyst,  but 

also  because  it  provides  reasonable  models  for  catalyst-support  interactions 

and  the  effects  of  sulfur  "poisoning."  The  biocatalytical ly  active  iron 

cluster  embedded  in  a cage  formed  by  the  neighboring  sulfur  atoms  and  the 

surrounding  organic  ligands  (for  the  detailed  structure,  see  Fig.  8)  is 

somewhat  analogous  to  the  embedding  of  a catalytically  active  iron  aggregate 

in  a supporting  refractory  material,  typically  silica  or  alumina  in  the  case 

15 

of  a commercial  catalyst.  Ironically,  sulfur  is  a strong  poison  for 

64 

commercial  iron  catalysts,  and  the  different  roles  it  play  electronically 

in  the  protein  and  in  the  catalyst  are  of  great  interest. 

The  structure  u of  the  4-iron-sulfur  protein  analogue  [Fe^SJfSCH^] 

is  shown  in  Fig.  8.  It  consists  of  a simple  cube  with  four  Fe  atoms  and 

four  "inorganic"  S*  atoms  placed  on  alternate  vertices.  Along  each  of  the 

diagonals  containing  the  Fe  atoms,  there  is  one  "organic"  S atom  terminated 

by  a methyl  group  CH^.  The  ideal  overall  point  symmetry  is  T^  and  the 

average  interatomic  distances  are  also  shown  in  Fig.  8.  However,  detailed 

X-ray  diffraction  studies^0  indicate  that  the  cube  is  slightly  distorted 

into  a lower  °2d  symmetry,  that  the  Fe-S-C  angle  is  not  exactly  180°,  and 

that  the  symmetric  interatomic  distances  are  not  exactly  equal.  All  these 

deviations  are  treated  as  perturbations  to  the  more  symmetric  model. 

2-  35 

The  orbital  energies  of  [Fe^S^SCH^)^]  , calculated  by  Yang  using 

the  SCF-Xa-SW  method  and  labeled  according  to  the  irreducible  representations 
of  the  T^  point  group,  are  displayed  in  Fig.  9.  The  highest  occupied 
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orbital  Is  Indicated  by  an  arrow,  and  the  SCF-Xa  energy  levels  for  the  free 
Iron  and  sulfur  atoms  are  shown  for  comparison.  The  energy  scale  is  such 
that  the  top  of  the  manifold  of  predominantly  sulfur-like  (S*-Fe)  orbitals 
coincides  in  energy  with  the  sulfur  3p  atomic  orbital.  The  following  points 
emerge  from  this  calculation: 

(1)  The  number  and  character  of  occupied  orbitals  suggest  that  each  Fe 

2 5+ 

atom  has  the  effective  fractional  valence  state  Fe  ' ; discrete 

Integral  valence  states  are  not  spectroscopically  detectable  in 
[Fe4S$(SCH3)4]2'.65 

(2)  The  electronic  inequivalence  of  S*  and  S atoms  is  emphasized  by  the 
findings  that  the  charge  localized  in  sulfur  atomic  regions  averages 
54.3  and  36.7%  in  the  S*-Fe  and  Fe-S  orbitals  respectively,  a behavior 
attributable  to  the  influence  of  the  electropositive  metal  in  the 
structural  fragments  S*-Fe3  and  Fe-S-C  (cf.  Fig.  8).  Kramer  et  al . 
measured  the  binding  energies  of  the  core  electrons  of  reduced  HP 
protein  by  X-ray  photoelectron  spectroscopy  and  found  two  sulfur  2p 
peaks  at  162.9  eV  and  161.5  eV,  which  they  assigned  to  the  S and  S* 
atoms  respectively.  However,  the  SCF-Xa  results  suggest  that  a 
reversed  assignment  is  more  appropriate,  since  the  local  potential 
wells  are  deeper  for  the  S*  atoms  than  the  S atoms  because  of  stronger 
electropositive  metal  influence,  and  consequently  greater  core-electron 
binding  energies.  It  is  clear  from  Fig.  9 that,  in  general,  the 
binding  energies  of  the  S*  atoms  are  greater  than  those  of  the  S atoms 
(e.g. , the  sulfur  3s  levels). 

(3)  The  band  of  levels  between  -0.35  and  -0.15  Ry  corresponds  to  orbitals 
("Fe-Fe'')  predominantly  Fe  3d- 1 i ke  in  character  and  orbitals  ("Fe-S") 
with  almost  equal  amounts  of  Fc-3d  and  S-3p  character.  However,  the 
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Fe-fe  levels  are  not  purely  metallic  in  nature,  as  compared  with  the 
orbitals  of  an  isolated  tetrahedral  Fe^  cluster  (cf.  Fig.  11)  having 
the  same  Fe-Fe  interatomic  distance.  The  Fe-Fe  orbitals  in  [Fe^S^( SCH^ 
have  S 3p-like  components  ranging  from  only  3*  S and  1%  S*  in  the  most 
metallic  occupied  bonding  orbital  (3e)  to  as  much  as  25%  S and  12%  S* 

In  the  least  metallic  unoccupied  antibonding  orbital  (5e).  Therefore, 
the  Fe-S  Interactions  are  quite  covalent. 

(4)  The  highest  occupied  level  (K^)  in  [Fe^S^SCH^)^]  is  predominantly 
tetrametal  antibonding  in  character,  with  S and  S*  contributions  to  the 
orbital  charge  amounting  to  18  and  7%  respectively.  A contour  map  of 
the  lOt^  orbital  wavefunction  is  shown  in  Fig.  10(a),  plotted  in  a 
Fe2SJ)  plane  (cube  face).  This  map  reveals  the  3d- 1 i ke  lobes  centered 
on  the  Fe  nuclei,  polarized  by  in-plane  and  out-of-plane  covalent 
interactions  with  S*  and  S 3p- 1 i ke  components,  respectively.  Also 
shown  in  Fig.  10(b)  and  10(c)  are  the  heavily  tetrametal  8t2  and  3e 

orbitals  which  illustrate  the  net  Fe-Fe  bonding. 

2-57 

Optical  absorption  spectra  for  [Fe^fSCH^]  are  somewhat  red- 
shifted  but  otherwise  are  similar  to  those  of  the  proteins  (HPred  and 
Fdox)  and  are  dominated  by  intense  features  near  295  nm  (4.20  eV)  and 
418  nm  (2.97  eV),  additional  incompletely  resolved  bands  at  intermediate 
energies,  and  a low  intensity  shoulder  at  'v  650  nm  (1.91  eV).  According  to 
the  SCF-Xa  calculation  (Fig.  9),  the  intense  peaks  and  fine  structure 
between  these  peaks  can  be  assigned  as  "charge-transfer"  excitations  from 
the  band  of  predominantly  sulfur-like  levels  lying  between  -0.48  and  -0.58 
Ry  to  the  partially  unoccupied  tetrametal  level  1 0t2*  The  top  ( 6t£ ) and 
bottom  (4aj)  of  the  sulfur  band  are  separated  from  the  lOtj  level  by  the 
energies  2.98  and  4.32  eV,  which  are  remarkably  close  to  the  energies  2.97 
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and  4.20  eV,  of  the  principal  absorption  peaks.  (Transition-state  calculations 
for  these  excitations  have  also  been  performed  and  the  relaxation  effects  are 
small.)  The  low  Intensity  shoulder  at  n,  1.91  eV  is  assigned  to  "d  d" 
electronic  transitions  between  the  occupied  and  unoccupied  tetrametal  ("Fe-Fe") 
levels. 

o_ 

The  measured  structure  of  the  analogue  [Fe^S^SCH^]  exhibits  D>d 

symnetry6^  and  those  of  the  biological  counterparts  HPrgd  and  FdQX  have  D d 
67  6R 

or  lower  symmetry,  ’ ° a property  which  is  explicable  in  terms  of  the  Oahn- 

2 

Teller  effect.  From  the  calculated  electronic  structure  of  the  [Fe^S^(SCH^)^] 

cluster  (Fig.  9),  the  degeneracy  of  the  partially  occupied  lOt^  orbital 

suggests  the  likelihood  of  a Jahn-Teller  distortion  from  the  perfectly  cubic 

(Tj)  geometry  to  one  of  lower  symmetry.  SCF-Xct-SW  calculations  for  both  the 

Td  and  °2d  structures  indicate  a lower  total  energy  for  the  distorted  cube 

and  a splitting  of  the  lOtg  level  into  e and  b£  levels,  as  shown  in  Fig.  9(b). 

This  splitting  amounts  of  approximately  0.005  Ry  (0.068  eV).  The  resulting 

"closed-shell"  (zero-net-spin)  electron  configuration  is  consistent  with  the 

observed  low-temperature  diamagnetism  of  the  analogue  and  HP  ^ The 

paramagnetism  observed  at  higher  temperatures6^ can  be  explained  in  the 

present  model  by  the  thermally  induced  excitation  of  electrons  across  the 

small  c ■+  bg  energy  gap  and  the  population  of  b2 1 states.  Lastly,  this 

model  suggests  that  the  oxidized  species  [Fe^S^SCH^]1 ~ and  its  biological 

counterpart  HPQx  would  have  the  (lOtg)  configuration  (a  half-filled  shell) 

and  hence,  in  the  spin-unrestricted  orbital  description,  would  not  be  subject 

to  Jahn-Teller  distortions.  Strong  paramagnetism  (S  = ^)  should  also  be 

observed  at  low  temperatures.  No  synthetic  analogue  of  this  oxidation  state 

has  been  produced  yet.  However,  it  has  been  observed  that  in  the  biological 

67 

state,  the  HPQx  site  is  marginally  more  symmetric  than  the  HP  d site, 
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and  that  the  cluster  volume  contraction  in  the  oxidation  process  HP^  -»  HPQx 

Is,  as  originally  suggested,^  partially  accountable  in  terms  of  electron 

removal  from  an  Fe-Fe  antibonding  orbital  (the  10t2  orbital  or  its  components 

in  a lower  symmetry).  It  has  also  been  reported  that  the  4-Fe-S  site  in 

the  Mo-Fe  protein  isolated  from  "Azatobacter  vinelandii"  has  a paramagnetic 

3 

ground  state  corresponding  to  S = j-. 

In  order  to  analyze  the  effects  of  the  "supporting"  sulfur-organic 
environment  on  the  electronic  structure  of  the  4-Fe  aggregate,  in  Fig.  11 
we  compare  the  SCF-Xa  energy  eigenvalues  of  the  predominantly  tetrametal 
(Fe-Fe)  orbitals  of  the  [Fe^S^SCH^)^]^’  cluster  (extracted  from  Fig.  9) 
with  the  spin-restricted  orbital  energies  of  a tetrahedral  Fe^  cluster  in 
which  the  Fe-Fe  internuclear  distances  have  been  constrained  to  the  values 
in  the  4-Fe-S  cluster  (cf.  Fig.  8).  The  following  conclusions  can  be  derived 
from  this  comparison: 

(1)  Tetrametal  bonding  in  the  elemental  Fe^  cluster  involves  both  d orbitals 
and  s,p  orbitals,  in  analogy  to  the  overlap  of  the  d band  by  the  s,p 
"conduction"  band  in  bulk  iron,  whereas  in  the  protein  analogue  the 
tetrametal  bonding  is  due  entirely  to  the  d orbitals. 

(2)  The  high  degeneracy  around  the  Fermi  level  of  the  Fe^  cluster  is 
removed  by  "turning  on"  the  sulfur-organic  environment  in  the  protein 
analogue. 

(3)  More  antibonding  tetrametal  d orbitals  are  occupied  relative  to  bonding 
orbitals  in  the  Fe^  cluster  than  in  the  protein  analogue,  Indicating 
that  the  tetrairon  d-orbital  bonding  is  actually  enhanced  by  the  sulfur- 
organic  environment. 

The  total  tetrametal  d-band  width  in  the  protein  analogue  is  significantly 
larger  than  that  of  the  elemental  Fe^  cluster,  consistent  with  the 


(4) 
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stabilization  of  the  bonding  by  the  sulfur-organic  environment  implied 
In  (3). 

An  even  more  dramatic  comparison  can  be  made  between  the  tetrairon 

d-orbltal  manifold  in  the  iron-sulfur  orotein  analogue  (Fig.  11)  and  the 

minority-spin  ( + ) d-orbltal  manifold  in  the  Fe^  cluster  analogue  of  a small 

Iron  particle  (Fig.  7),  In  relation  to  the  LUMO  of  the  nitrogen  molecule. 

This  provides  a common  basis  for  understanding  the  biocatalytic  mechanism 

of  nitrogen  fixation  (reduction  of  dinitrogen  to  ammonia)  by  the  enzyme 
33 

nitrogenase  and  the  heterogeneous  catalytic  mechanism  of  ammonia  synthesis 
29  30 

(the  Haber  process)  * by  a commercial  iron  catalyst.  As  argued  in 

Section  IV,  the  local  spin  polarization  of  pure  iron  aggregates  contributes 

to  the  reduction  in  activation  energy  (as  compared  with  other  metals)  for  the 

dissociation  of  Np,  the  rate-limiting  step  of  ammonia  synthesis,  by  reducing 

the  energy  (and  electronegativity)  difference  for  electron  flow  between  the 

catalyst  HOMO  and  Np  LUMO.  In  biological  nitrogen  fixation,  it  is  speculated 

that  the  active  center  of  the  4-Fe-S  protein  in  nitrogenase  reduces  dinitrogen 

33 

which  Is  bound  to  an  Mo-Fe  protein.  In  analogy  to  the  suggested  mechanism 
for  N2  dissociation  by  an  iron  catalyst  (cf.  Figs.  7 and  11),  we  see  that 
this  reduction  can  occur  by  electron  transfer  from  the  HOMO  of  the  tetrairon 
d band  of  the  4-Fe-S  cluster  to  the  LUMO  of  Np  (neglecting  the  perturbation 
of  this  LUMO  by  the  Mo-Fe  protein).  Furthermore,  the  tetrairon  HOMO  is 
significantly  closer  In  energy  and  electronegativity  to  the  Np  LUMO  than  in 
the  case  of  Fe^  and  larger  iron  aggregates,  suggesting  a lower  activation 
energy  for  electron  transfer.  This  seems  to  be  consistent  with  the  empirical 
fact  that  biological  nitrogen  fixation  occurs  efficiently  at  ordinary  tempera 
tures,  whereas  commercial  ammonia  synthesis  is  a high -temperature  process. 

This  theoretical  model  is  further  supported  by  the  observation  that  the 


Johnson  - 33 


53 

active  sites  in  nitrogenase  are  extremely  sensitive  to  oxygen  exposure, 
suggesting  the  mechanism  for  O2  bond  cleavage  and  oxidation  described  for 
metallic  iron  in  Section  IV.  Thus,  while  local  spin  polarization  gives 
metallic  iron  an  electron  configuration  which  promotes  the  catalytic  dissocia- 
tion of  N2  and  synthesis  of  ammonia,  it  is  the  supporting  sulfur-organic 
environment  which  stabilizes  the  tetrairon  cluster  in  an  electron  configura- 
tion highly  favorable  for  dinitrogen  reduction. 

The  final  issue  to  be  resolved  is  the  apparent  contradiction  between 

the  observation  that  sulfur  effectively  acts  as  a promoter  of  biocatalytic 

activity  in  iron-sulfur  proteins,  whereas  it  is  well  known  to  be  a severe 

64 

poison  of  iron  catalysts  in  anmonia  synthesis.  An  examination  of  the 

o 

electronic  charge  distribution  in  the  [Fe^S^SCH.^]  analogue,  as  calculated 
by  the  SCF-Xa  method,  suggests  that  the  Fe(3d)-S(3p)  and  Fe(3d)-S*(3p)  inter- 
actions in  the  corresponding  protein  are  primarily  covalent.  This  covalency 
is  fairly  apparent  in  the  orbital  contour  maps  shown  in  Fig.  10  and  is 
confirmed  experimentally  by  NMR  contact  shifts  measured  in  the  protein 
analogue. Consequently,  the  tetrametal  d orbitals  are  pushed  closer  in 
energy  and  electronegativity  to  the  Ng  LUMO,  as  compared  with  those  of  an 
elemental  iron  cluster  (cf.  Figs.  11  and  7),  because  of  the  antibonding 
components  of  the  iron-sulfur  covalency.  In  contrast,  the  interaction  of 
sulfur  with  metallic  iron,  as  represented  by  the  chemisorption  of  sulfur  on 
an  iron  surface,  is  more  Ionic,  withdrawing  electrons  from  the  metal  to  sulfur, 
thereby  increasing  the  work  function  of  the  metal  as  well  as  the  difference 
in  energy  and  electronegativity  between  the  iron  HOMO  and  N2  LUMO.  Thus 
coordinative  saturation  of  active  iron  surface  sites  by  sulfur  atoms, 
combined  with  the  changes  in  effective  iron  orbital  electronegativity,  leads 
to  a poisoning  of  iron  catalysts  for  Ng  dissociation  and  ammonia  synthesis. 


Johnson  - J4 


Despite  the  wide-ranging  similarities  of  the  measured  physical  properties 

p. 

of  analogues  such  as  [Fe^S|(SCH3)4]  to  those  of  the  active  centers  of  the 
4-Fe-S  proteins ^7,60-62,65  prote,n  environment  seems  to  be  necessary 
for  catalytic  activity.  Attempts  to  duplicate  the  full  process  of  nitrogen 
fixation  in  the  laboratory,  using  simpler  synthetic  analogues,  have  not  yet 
been  successful,  although  research  in  this  area  is  important  and  continues."^ 

It  is  not  yet  clear  to  what  degree  the  extended  organic  environment  of  the 
active  4-Fe-S  centers  in  the  enzyme  influences  the  local  electronic  structure 
and  chemical  properties  of  these  sites,  although  SCF-Xct-SW  studies^1  of  the 
substitution  of  other  simple  organic  ligands  (including  hydrogen)  for  the 
CH3  moiety  In  [Fe^SJfSCH^]  " do  not  qualitatively  alter  the  main  features 
of  the  electronic  structure  shown  in  Fig.  9 and  the  conclusions  based  thereon. 

VI.  CONCLUSIONS 

In  this  paper,  we  have  attempted  to  elucidate  a number  of  important 
analogies  among  the  electronic  structures  of  catalytically  active  sites  on 
transition-metal  surfaces,  in  isolated  transition-metal  coordination  complexes, 
in  supported  metal  catalysts,  and  in  metalloenzymes,  using  the  results  of 
first-principles  SCF-Xot  molecular-orbital  calculations  in  conjunction  with 
the  concept  of  spin-orbital  electronegativity.  Although  an  effort  has  been 
made  to  choose  realistic  illustrative  examples  for  which  there  are  comparable 
experimental  data,  some  of  the  conclusions  presented  herein  are  speculative 
and  hopefully  will  stimulate  further  experimental  and  theoretical  work. 

It  is  hoped  that  such  fundamental  interdisciplinary  studies  will  not  only 
lead  to  a better  understanding  of  the  chemical  properties  of  surfaces  but 
will  also  serve  ultimately  as  a guide  for  optimizing  catalytic  activity  and 
selectivity  through  the  systematic  refinement  of  existing  catalytic  materials 
or  through  the  molecular  design  of  entirely  new  catalysts. 
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FIGURE  CAPTIONS 

SCF-Xa  orbital  energies  for  coordinatively  unsaturated  transition- 
metal  complexes  representing  low-coordination  transition-metal 
sites  and  dissociative  hydrogen  chemisorption  thereon.  The 
highest  occupied  orbital  is  indicated  by  the  "Fermi  level"  c^. 
Contour  maps  of  the  principal  bonding  and  antibonding  molecular- 
orbital  wavefunctions  corresponding  to  the  orbital  energies  of 
the  LgMHg  complex  shown  in  Fig.  1. 

Possible  reaction  path  for  the  hydrogenation  of  acetylene  at  a 
coordinatively  unsaturated  transition-metal  site. 

Spin-polarized  SCF-Xa  orbital  energies  of  a 9-atom  bcc  iron 
cluster. 

Spin-polarized  SCF-Xa  orbital  energies  of  a 15-atom  bcc  iron 
cluster. 

Contour  maps  of  principal  bonding  and  antibonding  spin  orbitals 
of  bcc  Fe^5  cluster.  Solid  curves  indicate  positive  parts  of 
wavefunction;  dashed  curves  indicate  negative  parts  of  wavefunc- 
tion.  e^  orbitals  mapped  in  (100)  plane  containing  central 
atom  and  second  nearest  neighbors;  t^g  orbitals  mapped  in  (110) 
plane  containing  central  atom,  nearest-neighbor  atoms,  and 
second  nearest  neighbors;  t^g  orbital  mapped  in  plane  of  cube 
face. 

Comparison  of  the  SCF-Xa  orbital  energies  of  the  Fe^  and  Pt^ 
clusters  with  those  of  the  N^,  CO,  and  0^  molecules.  Also  shown 
for  comparison  are  the  SCF-Xa  atomic-orbital  energies  of  N,  C, 
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Figure  8.  Structure  of  [Fe^S^SCH^]  , a synthetic  analogue  of  the  active 
center  of  oxidized  ferredoxin  and  reduced  high-potential  four- 
iron-sulfur  proteins. 

Figure  9.  (a)  SCF-Xa  orbital  energies  of  [Fe^S|(SCH3)^]  in  symmetry 

and  (b)  splitting  of  highest  occupied  orbital  (10t2)  in  D2(J 
symmetry.  SCF-Xa  atomic-orbital  energies  for  Fe  and  S are  also 
Included  for  comparison. 

Figure  10.  Contour  maps,  plotted  in  Fe2S£  plane  (cube  face),  for  the  (a) 

10t2,  (b)  8t2,  and  (c)  3e  orbital  wavefunctions  of  [Fe^S^SCH^] 

Figure  11.  Comparison  of  SCF-Xa  orbital  energies  for  an  elemental  tetrahedral 
Fe4  cluster  with  the  corresponding  tetrametal  orbital  energies 
for  the  protein  analogue  [Fe^S^SCH^]  . Also  shown  for 
comparison  are  the  SCF-Xa  orbital  energies  of  the  N2  molecule. 
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